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Abstract 

The fundamentals of a quantum heat engine are derived from first principles. The study is based 
on the equation of motion of a minimum set of operators which is then used to define the state 
of the system. The relation between the quantum framework and thermodynamical observables 
is examined. A four stroke heat engine model with a coupled two- level-system as a working fluid 
is used to explore the fundamental relations. In the model used, the internal Hamiltonian does 
not commute with the external control field which defines the two adiabatic branches. Heat is 
transferred to the working fluid by coupling to hot and cold reservoirs under constant field values. 
Explicit quantum equation of motion for the relevant observables are derived on all branches. The 
dynamics on the heat transfer constant field branches is solved in closed form. On the adiabats, 
a general numerical solution is used and compared with a particular analytic solution. These 
solutions are combined to construct the cycle of operation. The engine is then analyzed in terms 
of frequency-entropy and entropy-temperature graphs. The irreversible nature of the engine is the 
result of finite heat transfer rates and friction-like behavior due to noncommut ability of the internal 
and external Hamiltonian. 
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I. INTRODUCTION 



Analysis of heat engine models has been a major part of thermodynamic development. 
For example Carnot's engine preceded the concepts of energy and entropy 1]. Szilard and 
Brillouin constructed a model engine which enabled them to resolve the paradox raised by 

nn 

Maxwell's demon |2J, |3(. The subsequent insight enabled the unification of negative entropy 
with information. In the same tradition, the present paper studies a heat engine model with 
a quantum working fluid for the purpose of tracing the microscopic origin of friction. The 
function of a quantum heat engine as well as its classical counterpart is to transform heat 
into useful work. In such engines, the work is extracted by an external field exploiting the 
spontaneous flow of heat from a hot to a cold reservoir. The present model performs this 
task by a four stroke cycle of operation. All four branches of the cycle can be described by 
quantum equations of motion. The thermodynamical consequences can therefore be derived 
from first principles. 

The present paper lays the foundation for a comprehensive analysis of a discrete model 
of a quantum heat engine. A brief outline which has been published emphasized the engines 
optimal performance characteristics j^J. It was shown that the engines power output vs. 
cycle time mimics very closely a classical heat engine subject to friction. The source of the 
apparent friction was traced back to a quantum phenomena: the noncommutability of the 
external control field Hamiltonian and the internal Hamiltonian of the working medium. 

The fundamental issue involved require a detailed and careful study. The approach 
followed is to derive the thermodynamical concepts from quantum principles. The connecting 
bridges are the quantum thermodynamical observables. Following the tradition of Gibbs a 
minimum set of observables is sought which are sufficient to characterize the performance 
of the engine. When the working fluid is in thermal equilibrium, the energy observable is 
sufficient to completely describe the state of the system and therefore all other observables. 
During the cycle of operation the working fluid is in a non-equilibrium state. In frictionless 
engines, where the internal Hamiltonian commutes with the external control field, the energy 
observable is still sufficient to characterize the engine's cycle 0, f|. In the general case 
additional variables have to be added. For example in the current model, a set of three 
quantum thermodynamic observables is sufficient to characterize the performance. With 
only two additional variables the state of the working fluid can be characterized also . A 
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knowledge of the state is necessary in order to evaluate the entropy and the dynamical 
temperature. These variables are crucial in establishing a thermodynamic perspective. 



re current investigation is in line with previous studies of quantum heat engines 



SHH0, 0,00,0 



17. 
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181 ] . All the studies of first principle quantum models 



have conformed to the laws o 



resembling turbines 



thermodynamics. These models have been either continuous 



17j . or discrete as in the present model 
Surprisingly the performance characteristics of the models were in close resemblance to 
their realistic counterparts. Real heat engines operate far from the reversible conditions, 
where the maximum power i s restricted due to finite heat transfer 0] , internal friction and 



heat leaks 120 
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26| . Analysis of the quantum models of heat engines, based 
on a first principle dynamical theory, enable to pinpoint the fundamental origins of finite 
heat transfer, internal friction and heat leaks. 

Studies of quantum continuous heat engine models have revealed most of the known 
characteristics of real engines. In accordance with finite time thermodynamics the power 
always exhibits a definite maximum [2^, and the performance has been limited by heat leaks 



12 1 . Finally indications of restrictions due to friction like phenomena have been indicated 
12j. The difficulty with the analysis is that it is very hard to separate the individual 



contributions in the case of a continuous operating engine. 

To facilitate the interpretation a four stroke discrete engine has been chosen for analysis. 
The cycle of operation is controlled by the segments of time that the engine is in contact 
with a hot and cold bath and by the time interval required to vary the external field. To 
simplify the analysis the time segments where the working fluid is in contact with the heat 
baths are carried out at constant external field. Such a cycle of operation resembles the 
Otto cycle which is composed of two isochores where heat is transfered and two adiabats 
where work is done. This simplification allows to obtain the values of the thermodynamical 
observables during the cycle of operation from first principles in closed form. 



II. QUANTUM THERMODYNAMICAL OBSERVABLES AND THEIR DYNAM- 
ICS 

The quantum thermodynamical observables constitute a set of variables which are suf- 
ficient to completely describe the heat engine performance characteristics as well as the 
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entropy and temperature changes of its working medium. The analysis of the performance 
requires a quantum dynamical description of the changes in the thermodynamical observ- 
ables during the engine's cycle of operation. The thermodynamical observables are associ- 
ated with the expectation values of operators of the working medium. Using the formalism 
of von Neumann, an expectation of an observable (A) is defined by the scalar product be- 
tween the operator A representing the observable and the density operator p representing 
the state of the working medium: 

(A) = (A • p) = Tr{At p} . (1) 

The dynamics of the working medium is subject to external change of variables as well 
as heat transport from the hot and cold reservoirs. The dynamics is then described within 



the formulation of quantum open systems 



28, 



where the dynamics is generated by the 



Liouville super operator C either as an equation of motion for the state p (Schrodinger 
picture): 

P = C(p) , (2) 
or as an equation of motion for the operator (Heisenberg picture): 

A = £*(A) + °— . (3) 



The second part of the r.h.s. appears since 



;he operator A can be explicitly time dependent. 



Significant simplification is obtained |27( when: 
• a) The operators of interest form an orthogonal set Bi i.e. 

% ■ Bj) = 5 ;j , (4) 
where Bo = I is the identity operator. 

b) The set is closed to the operation of C* . 

i = £*(B0 = ^jBj , (5) 
j 

where l\ are scalar coefficients composing the matrix L. 

c) The equilibrium density operator is a linear combination of the set: 

P e9 = ^ i + E^B k , (6) 



N 

k 



where N is the dimension of the Hilbert space and are the equilibrium expectation 
values of the the operators, (B^ 3 ) . 

The operator property of Eq. (J5J) allows a direct solution to the Heisenberg equation of 
motion (jSJ) by diagonalizing the L matrix, relating observables (Bk) at time t to observables 
at time t + At that is h(t + At) = U(At)h(t) where U = e LAt and b is a vector composed 
from the expectation values of B k ( for an example Cf. (J35j) ). 

The time dependent expectation values b(t) and Eq. (0) can be employed to reconstruct 
the density operator: 

PR = ^ i + E^B k , (7) 

k 

where the expansion coefficients become b^ = (Bk). Although the set Bk is not necessarily 
complete, equation ((7J) will still be used as a reconstructing method for the density operator. 
This reconstructed state p R reproduces all observations which are constructed from linear 
combinations of the set of operators B^. 

The Liouville operator Eq. (0),© for an open quantum system can be partitioned into 
a unitary part C H and a dissipative part C D [28]: 

C = £ H + C D . (8) 

The unitary part is generated by the Hamiltonian: H: 

Z^(A)=i[H,A] . (9) 



The condition for a set of operators to be closed under C* H have been well studied |29J. If 
the Hamiltonian can be decomposed to: 



H = $>jBj , (10) 



and the set B k forms a Lie algebra |30j, |3l| i.e. [B i; Bj] = C;jB k (the coefficients are 
the structure factors of the Lie algebra), then the set is closed under C*„ . 



For the dissipative Liouville operator Co, Lindblad's form is used 
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C* D {k) = (FjAFj - ^(FjFjA + AFjF})) , (11) 

j 

where Fj are operators from the Hilbert space of the system. The conditions for which the 
set Bj is closed to C* D have not been well established. Nevertheless in the present studied 
example such a set has been found. 



A. Energy balance 



The energy balance of the working medium is followed by the changes in time to the 
expectation value of the Hamiltonian operator. For a working medium composed of a gas 
of interacting particles the Hamiltonian is described as: 

H = H ext + H int . (12) 

H ex t — is the sum of single particle Hamiltonians, where u = u(t) is the time de- 

pendent external field. It therefore constitutes the external control of the engine's operation 
cycle. Hi n t represents the uncontrolled inter-particle interaction part. 

The existence of the interaction term in the Hamiltonian means that the external field 
only partly controls the energy of the system. One can distinguish two cases, the first is 
when the two parts of the Hamiltonian H ext and Hi nt commute. The other case occurs 
when [H ext , Hi n t] 7^ leads to [Hi nt (t), Hi nt (t')] 7^ 0, causing important restrictions on the 
cycle of operation (Cf. section |VHj). 

Since the energy is E = (H), the energy balance becomes Cf. Eq. (J3J): 

f - < £ -<«» + <i > ■ <"> 

Eq. (fT3J) is composed of the change in time due to the explicit time dependence of the 
Hamiltonian (Cf. Eq. © interpreted as the thermodynamic power: 

V = ^<Hi> , (14) 

i 

where (Hi) is the expectation value of the single particle Hamiltonian. The accumulated 
work on an engines trajectory W = J Vdt. 

The heat flow represents the change in energy due to dissipation: 

Q = (£b(H)) = (£* D (H ext + Hint)) , (15) 

(note £*(H) = £* D (H) since C* H {H) = 0). Egs. (E L dH) and (JSD leads to the time 
derivative of the first law of thermodynamics 0, [lg, 3| : 

§ - v + a . (i6) 
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B. Entropy balance 



Assuming the bath is large the entropy production due to heat transfer from the system 
to the bath becomes: 

VS = I , (17) 

where T is the bath temperature. 

Adopting the supposition that entropy is a measure of the dispersion of the measurement 
of an observable (A), we can label the entropy of the working medium according to the 
measurement applied i.e. S^. The probability of obtaining a particular ith measurement 
outcome is: Pi = tr{Pip} where Pj = are the projections of the i th eigenvalue of the 
operator A. The entropy associated with the measurement of A becomes: 

$a = -J^logPi , (18) 

i 

The probabilities in Eq. (|18|) can be obtained from the diagonal elements of the density 
operator p in the eigen-representation of A. The entropy of the operator A that leads 
to minimum dispersion (|18jh defines an invariant of the system termed the Von Neumann 



entropy 



35j: 



Svn = -tr{p\ogp} , (19) 



<?a — $vn f° r ah A. The analysis of the energy entropy Se = S^ of the working fluid 
during the cycle of operation is a source of insight into the dynamics. It has the property: 
Se > <Svn with equality when the p is diagonal in the energy representation which is true 
in thermal equilibrium. Then: 

e -/3H 

P eq = — , (20) 

with (3 = 1/kfjT and Z = tr{e _/3H }, The systems temperature has thus become identical 
with the bath temperature. When the working medium is not in thermal equilibrium, a 
dynamical temperature of the working medium is defined by j^]: 

(M) 

Tdyn = fdS E \ ' (21) 
V dt ) 

and will be used to define the internal temperature of the working fluid (Cf. Section IV)) . 
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III. THE QUANTUM MODEL 



The following quantum model demonstrates a discrete heat engine with a cycle of opera- 
tion defined by an external control on the Hamiltonian and by the time duration where the 
working medium is in contact with the hot and cold bath. The model studied is a particular 
realization of the general framework of section |H] First the generators of the motion Lh and 
Co are derived leading to equations of motion. These equations of motion are then solved 
for each of the branches thus constructing the operating cycle. 

A. The equations of motion 

The generators of the equations of motion are the Hamiltonian for the unitary evolution 
and Cn for the dissipative part (Cf. Eq. JHJ)). 

1. The Hamiltonian 

The single particle Hamiltonian is chosen to be proportional to the polarization of a two- 
level-system (TLS): & 3 Z , which can be realized as an ensemble of spins in an external time 
dependent magnetic field. The operators & z ,& x , & y are the Pauli matrices. For this system, 
the external Hamiltonian, Eq. (fT^j) becomes: 

H ext = 2- 3 / 2 ^)(^®i 2 + i 1 ®<rf) , (22) 

and the external control field uo(t) is chosen to be in the z direction. The uncontrolled inter- 
action Hamiltonian is chosen to be restricted to coupling of pairs of spin atoms. Therefore 
the working fluid consists of noninteracting pairs of TLS's. For simplicity, a single pair can 
be considered. The thermodynamics of M pairs then follows by introducing a trivial scale 
factor. Accordingly the uncontrolled part is: 

Hint = 2-^j(d- l x ®6-l-crl®crl .) (23) 

J scales the strength of the interaction. When J — > 0, the model represents a working 
medium with noninteracting atoms p. The interaction term, Eq. (J23)) . defines a correla- 
tion energy between the two spins in the x and y directions. As a result, the interaction 
Hamiltonian does not commute with the external Hamiltonian Eq. (|22jh which is chosen to 
be polarized in the z direction. 

8 



2. The operator algebra of the working medium 



The maximum size of the complete operator algebra of two coupled spin systems is 16. A 
minimum set of operators closed to C* is sought which is sufficient as the basis for describing 
the thermodynamical quantities. First, a Lie algebra which is closed to the unitary evolution 
part is to be determined. To generate this algebra the commutation relations between the 
operators composing the Hamiltonian are evaluated ( Cf. Eq. (jlOj)). Defining: 

/ 1 \ 


\0 -I J 

where the tensor product eigenstates of a z and &\ are used for the matrix representation, 
termed the "polarization representation". 



(24) 



The second operator B 2 is: 



B; 



2- 3/2 (*J 



y) 



1 

72 



/o l\ 


\ 1 0/ 



(25) 



The commutation relation: [Bi, B 2 ] = y2iBs leads to the definition of B 

/ - i\ 

B 3 = 2- 3 / 2 (*l + 



1 

7! 



0000 
0000 

\ i ) 



(26) 



The set of operators Bi, B2, B3 form a closed sub-algebra of the total Lie algebra of the 
combined system. The Hamiltonian expressed in terms of the operators B 1; B 2 , B 3 becomes: 

/ uj J \ 

1 

a/2 

\ J -uj ) 



H 



cjBi + JB 2 



(27) 
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TABLE I: Multiplication table of the commutation relations [X, Y] of the operators Bi between 
themselves and with the Hamiltonian. 



X\Y 


Bi 


B 2 


B 3 


Bi 





iV2B 3 


-iV2B 2 


B 2 


-iV2B 3 





iv^Bi 


B 3 




-i>/2Bi 





H 


-iv^JBs 




i\/2JBi - iV2ujB 2 



All the three operators are Hermitian, and orthogonal (Cf. Eq. (@J ). Table (I) summarizes 
the commutation relations of this set of operators. 

The commutation relations of the set of Bk operators define the SU(2) group and are iso- 
morphic to the angular momentum commutation relations by the transformation Bk — > Jk- 
Bi,B 2 ,B 3 can be identified as the generators of rotations around the z,x and y axes re- 
spectively. This representation allows to express the expectation values in a Cartesian three 
dimensional space ( See Fig. OJ) . 



3. The generators of the dissipative dynamics 

The dissipative part of the dynamics is responsible for the approach to thermal equi- 
librium when the working medium is in contact with the hot/cold baths. The choice of 

28| . The operators 



Lindblad's form in Eq. (fTT]) guarantees the positivity of the evolution 
Fj which lead to thermal equilibrium are constructed from the transition operators between 
the energy eigenstates. Diagonalizing the Hamiltonian (|T2"j) leads to the set of energy eigen- 
values and eigenstates: 

n n 

e 1 = - — , e 2 = 0, e 3 = 0, e 4 = -j= , (28) 



where Q = \/uj 2 + J 2 . The method of construction of Fj is based on identifying the operators 
with the raising and lowering operators in the energy frame. For example, Fi = \fk ||2)(1| 
or F 2 = \Jh ||1)(2|. The bath temperature enters through the detailed balance relation 
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HQ 



hi 

ki 



'75 



(29) 



The operators Fj constructed in the energy frame are then transformed into the polarization 
representation. The details are described in Appendix B. 
Substituting the Bi operators into £d, Eq. (JUJ), one gets: 



£d(Bi) 
£d(B 2 ) 
£d(B s ) 



-r(Bi + ^i) 



r(B 2 + 



J fcj-fcT 



>2 1 V2U r 

r(B 3 ) , 



(30) 



where T = k j + k T- 

From Eq. (J3*U|) the set of {B} operators and the identity operator I are invariant to the 
application of the dissipative operator Cd which leads to equilibration. 

The interaction of the working medium with the bath can also be elastic. These encoun- 
ters will scramble the phase conjugate to the energy of the system and are classified as pure 
dephasing (T 2 ) (Cf. Eq. (fTH|) . In Lindblad's formulation the dissipative generator of elastic 
encounters is described as: 



Z^e(A) 



-7[H, [H, A] 



(31) 



The elastic property is equivalent to ££> e (H) = 0. Moreover the set B ; which is closed to 



the commutation relation with H is also closed to £ 



To summarize the set B 1; B 2 , B 3 and I is closed under the operation of C* = C* H + C* D + 
C* De . Gathering together the various contributions leads to the explicit form of the equation 
of motion: 



d 
dt 



/ 



(B 2 ) 

V ^ J V 

or in vector form where bk 



r - 27 j 2 -27 ju V2J 

-T - 27a; 2 -\/2uj 

-V2J V2lu -r °~ n2 



2 7 fi 2 J 



( ( Bl ) N 

(B 2 > 
V (B 3 > j 



v4( fc i- fc t) (P2) 



(B k ): 



dt 



Bb 



(33) 



11 



B. Integrating the equations of motion 

The thermodynamical observables require the solution of the equations of motion on 
all branches of the engine. The field values uj are time independent on the isochores thus 
allowing a closed form solution, u changes with time on the adiabats therefore solving the 
equation of motion either requires a numerical solution or finding a particular solution based 
on an explicit time dependence of to. 



1. Solving the equations of motion on the isochores. 

On the isochores the coefficients in Eq. (|33|) are time independent. A solution is found 
by diagonalizing the B matrix leading to the eigenvalues: — Y — i\/2VL — 2^VL 2 , — Y 
and — T + i\^2il — 2 / yil 2 . The diagonalization enables to perform in closed form the 
exponentiation of e B At obtaining the propagator of the working medium operators Lf(At). 

' e -(r+iV2n+2jn 2 )At q q \ 

U(At) =11 e(- rAt ) 7?r J 



V 











~{T-iy/2tt+2~fn 2 )At 



where: 



K 



leading to the final result: 



( U/y/2Q co/n -u/V2~n\ 

-iu/V2Q J /ft iu/V2Q 
l/y/2 l/y/2 



U(At) = exp -(r + 2 7 fi 2 )At 



/ Xuj 2 +cJ 2 mJ(X-c) Js \ 

a, 2 n 2 n 

UJJ(X-C) XJ 2 +CLU 2 -LUS 



\ 



n 2 



n 2 

us 

n 



n 
c 



/ 



(34) 



(35) 



where X = exp(27^ 2 At), c = cos(v / 2^At) and s = sin( v / 2^At). The solution of Eq. (f32jl 
then becomes: 



b(t + At) = W(At)(b(t) -b e i) +b e( i , 



(36) 
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where the equilibrium values of the operators are calculated from the steady state solutions 
ofEq. fl33J): 

h(n/?/V2) = - 

i(np/V2) = -^^i^ (37) 

On the isochores the solution of Eq. (J35j) can be extended to the full duration Th/ C of 
propagation on the hot/cold branches. Therefore, At = 7"w c . 

There are cycles of operation where the external field u also varies when the working 



bf = 


( B l q ) = 


nz 




(B e 2 q ) = 


%/2J 

nz 


h eq - 
°3 — 


(B e 3 q ) = 


. 



medium is in contact with the hot or cold baths, for example the Carnot cycle For 
such cycles the equation of motion can be solved by decomposing these branches into small 
segments of duration At. Then Eq. (136(1 can be used as an approximate to the short time 
propagator. 



C. Propagation of the observables on the adiabats 

The equations of motion on the adiabats have explicit time dependence. To overcome this 
difficulty two approaches are followed. The first is based on decomposing the evolution to 
short time segments and using a short time approximation to solve the equations of motion. 
The second approach is based on finding a particular time dependence form of u(t) which 
allows an analytic solution. 



1. Short time approximation 



For the adiabatic branches the working medium is decoupled from the baths so that the 
time propagation is unitary. Eq. (|32|) thus simplifies to: 

/ 

d 



dt 







V2J \ ( O 



V 



/ 





-V2u(t) b 2 ■ (38) 

V2J y/2u(t) 

L(t)b. Since the matrix L(t) is time dependent the 
propagation is broken into short time segments At, reflecting the fact that [L(t), L(t')] ^ 0, 

3 At 



V h ) 



Or in the vector form: -jrb 

at 



b(t) 



N 

JJexp 



L{t')dt' b(0) 



0-1) At 



(39) 
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where NAt = t. Eq. ()38|) is solved by diagonalizing the matrix L for each time step assuming 
that during the period At uj(t) is constant. Under such conditions U a (t, At) becomes: ( the 
index a stands for adiabat) 



Ua(t,At) 



JL(t)M 



( 


UJ 2 +CJ 2 


wJ(l-c) 


Js 








Q 2 


n 






uj(l-c) 


J 2 +cw 2 








Q 2 


n 2 


n 






Js 




c 


J 


\ 


n 


n 



(40) 



which becomes the short time propagator for the adiabats from time t to t + At. 



2. An analytical solution on the adiabats 

The analytic solution for the propagator on the adiabats is based on the Lie group struc- 
ture of the {B} operators. The solution is based on the unitary evolution operator U(t) 
which for explicitly time dependent Hamiltonians is obtained from the Schrodinger equation: 

-ijV(t) = H(t)U(*), U(0)=i . (41) 
The propagated set of operators becomes: 

h{t) = V{t)h(<d)t\t) = U a (t)S(0) , (42) 
and is related to the super-evolution operator U a {t). Based on the group structure Wei and 



Norman, 37( constructed a solution to Eq. (j41j) for any operator H which can be written 
as a linear combination of the operators in the closed Lie algebra H(t) = X^j=i hj(t)Bi,, 
where the hi(t) are scalar functions of t, ( Cf. Eq. (jTUjl ). In such a case the unitary evolution 
operator U(t) can be represented in the product form: 

m 

U(t) = l[exp(a k (t)B k ) . (43) 

k=l 

The product form replaces the time dependent operator equation (|4(J|) with a set of scalar 
differential equations for the functions ctfc(t). As has been shown in 1111 A 2l three op- 
erators form a closed Lie Algebra. Writing the unitary evolution operator explicitly leads 
to: 

U(t) = exp^BOexp^B^exp^Bs) (44) 
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The a/2 factor is introduced for technical reasons. Based on the group structure [37] Eq. 
(j4*Tj) leads to the following set of differential equations has to be solved: 

/ \ T /Sin(ai) sin(a:2) s . n? T , % . V2Jsm(ai) 

«! = V2u (t) + V2J( — v ; . . v ; ! "2 = >/2Jcos(ai) ; d 3 = r^V^ ■ 45 

cos(a 2 ) cos(a 2 ) 

Using Eq. P2"j) the propagator U a {t) is evaluated explicitly in terms of the coefficients a: 



^ C2C3 — S3C1 + C3S2S1 C3S2C1 + S3S1 ^ 



W a (t) 



(46) 



C2S3 C3C1 + S3S2S1 S3S2C1 — C3S1 

y -s 2 C 2 Si C 2 Ci J 

where: S\ = sin^x), s 2 = sin(a 2 ), S3 = sin(a 3 ), Ci = cos(«i), c 2 = cos(a 2 ), c 3 = cos(a 3 ). 

The problem of obtaining a closed form solution for the propagator U a (t) has been trans- 
formed in to finding the solution of three coupled differential equations Eq. (J431 ) which 
depend on u(t). A general solution has not been found but by choosing a particular func- 
tional form for uj(t) a closed form solution has been obtained. 

3. The Explicit Solution for a 

To facilitate the solution of Eq. (J45|) . a particular form of uj{t) is chosen: 

u(t) =«L- j^M ™W . (47) 
V2 cos(a 2 ) 

Two auxiliary functions are defined, u{t) and v(t): 

u(t) = -J 2 t 2 + V2rJt- v(t) = r- V2Jt . (48) 

r is a constant which restricts the product Jt: { < r < 1; Jt < V2r}. In terms of u(t) 
and v(t), the solutions of Eq. (|43|) become: 

«i = arccos ( — = | . (49) 



«2 = arcsin ( — -(r\/l + 2w — v) ) . (50) 



a 3 = - - ln(2V4u 2 + 2u + 4u + 1) 
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Vl - r 2 r . fl. l-r 2 A . fl. 1 



r 



-{arcsin ( -fl 1 + arcsin [ -fl 1 } . (51) 

For t = 0, U = I, therefore ai(0) = 0, 02(0) = 0, 03(0) = which is consistent with Eq. 
AH,® and (JED). 

Introducing into Eqs. (J47j) the explicit functional forms of u(t) becomes: 

v / 2(l + 2u) v /u v / l + 2M(v / l + 2M + rt») 
At t = 0, is singular. Since the engine operates between two finite values of u a corre- 
sponding time segment is chosen which does not include the singularity at t — (Cf. Fig. |TJ). 
Using the group property of U a (t), i.e. U a (ti)U a (t2) = U a {t\ + t 2 ) the propagation is carried 
out by changing the origin of time, U a (t) = U~ l \to)U a {t + to) where to is either U for the 
compression adiabat or tf for the expansion adiabat. One should note, that U~ l {t) = U\{t) 
but due to the explicit time dependence U~ l {t) 7^ Ul(—t). 



IV. RECONSTRUCTION OF p R 



The reconstruction p is designed to describe the state of the working medium from its 
initial state to equilibrium. As was analyzed in the previous section, the set of operators 
B 1 ,B 2 ,B 3 ,I are sufficient to describe the energy changes during the cycle of operation of 
the engine. Is this set sufficient to reconstruct the density operator? 

In equilibrium p eq is diagonal in the energy representation, From the eigenvalues of the 
Hamiltonian Eq. ()28|) . p eq in the energy picture becomes: 



Z 









\ 











\ 



(53) 



where 
Z 



exp-({l{3/V2) + 2 + exp(0/3/V2) = -| + 2 + £i 

k J, K \ 



k Ik t 



(54) 
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FIG. 1: The external field w as a function of time on the adiabats corresponding to the function 
Eq. (|52|) for which an analytic solution exist. Indicated are the values of the initial and the final 
time and of the corresponding u which are used to construct the cycle of operation. Notice the 
singularity at t = 0. 

By inspection, the diagonal elements of the equilibrium density operator are seen to be 
denned by three independent variables. The energy expectation accounts for one variable. 
The expectation value of B3 has no diagonal elements in the energy representation therefore 
two additional operators are required to facilitate a reproduction of p R . 

/o o\ 



2 -3/2 



I 2 — I ] 



1 

7! 



0100 
00-10 
\o 0/ 



(55) 



and 



Br 



2*z ® 



( \ o\ 
0-100 
0-10 
\0 1/ 



(56) 



Since both B 4 and B 5 commute with the Hamiltonian, they undergo only dissipative dy- 
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namics but they are uninfluenced by the dephasing generated by C* D : 

*U = -TB 4 (57) 

with the solution: 

B 4 (t) = B 4 (0)exp(-rt) , (58) 
and (B^) = 0. The equation of motion of B 5 is: 

£ = -2rB 5 - >/2~(kl - fct)Bi- V2-(kl - fct)B 2 = 

-2rB 5 + 2r(Bi q )Bx + 2r(B2 q )B 2 . (59) 

At equilibrium, = 0, and then (Bg q ) = (B^ q ) 2 + (B^ 01 ) 2 , a result which can be 
verified by computing (Bg q ) = tr{p eg B5}. Eq. (|59|) is a linear first order inhomogeneous 
equation for B 5 depending on the time dependence of the closed set B 1 ,B 2 ,B 3 , Eq. 
Changing Eq. (|55|) to observables, Eq. (JTJ), and by integrating subject to the solutions of &i 
and 62 leads to: 

6s(t) = ? |(a;(6 1 (0)-&f) + J(6 2 (0)-6^))(a;6f + J6^)(e- r *-e- 2r *) + 



fco(fcic(e^ (r+ ^ )f ) + fc 2 se - (r+2 ^ 2) * - fc ie - 2r *) + (b 5 (0)-bt q )e~ 2rt + bt q , (60) 
where: 



2f( Jb eq - ujb eq ) 



(61) 



fi 2 ((r + 2 7 fi 2 ) 2 + 2^ 2 ) 

k x = (J(6 1 (0)-6f)-cu(6 2 (0)-6^))(r + 27fi 2 )-fi(6 3 (0)-^)(v / 2fi) 

and 

A; 2 = (J(&x(0) -6f)-o;(6 2 (0)-6f))(V2n) + fi(6 3 (0) - ^)(r + 2 7 fi 2 ) . 

Using the set of of the five orthogonal and normalized operators together with the identity 
operator the density operator p R is reconstructed. Representing p R in different basses 
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facilitates the calculation of the different entropies. p R in the polarization basis becomes: 



4 "T" ^ "T" 2 




I J>2_ _|_ „-^3_ 



4 ^ V2 
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_&2 „-^3_ \ 




1 bi i_ 65 

4 2 / 



(62) 



The off diagonal elements of p are the expectation values of the operators: 
B-t = -^(&2 =t 2B3) which represent the correlation between the individual spins. 
The density operator p R in the energy basis becomes: 



/ 1 E_ , 65 

4 fV2 ^ 2 










A _|_ »4 
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I 463 _ Jfei _|_ \ 






(63) 



n 1 4. E a. k / 

u 4 + + 2 / 

where = + J&2- In equilibrium, the off-diagonal elements vanish, and the matrix will 



be identical to Eq. ([53)1 . In non-equilibrium, the off-diagonal elements of p e determine the 
"phase" Cf. Sec. EEED 

To compute the Von-Neumann entropy p R is diagonalized leading to: 
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'4 y/2 + 2 
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(64) 



where D = ^Jh\ + h\ + bj. 



V. DYNAMICAL TEMPERATURE (T dyn ) ON THE BRANCHES 

Based on the definition of the dynamical temperature Td yn in Eq. (|2"T|) . and from Eq. 
(123): 



uh+uh + Jk uh-TE-^jki-k T) 

/;/ " Epf(i + iog(pf)) " -£*f(i + iog(pf)) * 



The four probabilities pf are the diagonal elements of the density operator in the energy 
representation p e , Eq. (JS3j). The derivatives of the probabilities are obtained from Eqs. 
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and dSSj) : 



. g _ -u;6i + TE | (fc | -fc t) , & 5 ._g _&5 
Pl " ^ + 2 + 2 ' P2 " 2 • 



. B _ &s . E uh-YE [kl-k]) b 5 



« = p ^^--^ + i • (66) 

where b 5 is obtained form Eq. (|55j): 65 = 2r(6^ 9 6i + &^&2 — &s)- 



A. Dynamical temperature on the isochores 

Evaluating the derivatives of the probabilities in Eqs. ()66|) and using the fact that on the 
isochores to =0, the dynamical temperature, Eq. (|65jl . becomes: 

(rs + ^(fc l -fc t)) 

T^ yn — —. ; r- . (67) 

(^log(piM) + (kl ~ 2 k]) log(piM) + |&5 l0g(piP4Mp 3 )J 

A consistency check is obtained by comparing T^ yn for J=0 with the internal temperature 
of a two- level-system. For J = 0: 

Tdyn = .^w( W^ ) ' (68) 

which leads to: 

" = " ' (69) 

which is the internal temperature for a noninteracting spin system with energy spacing uj/\/2 

0. 



B. Dynamical temperature on the adiabats 

On the adiabats 64 = and 65 = 0. From Eq. (|u3jl and (ju^jl the derivatives of the 
probabilities on the adiabats become: 

Pf = - ^; Pf = 0; pf n= ^ (TO) 

which leads to the dynamical temperature on the adiabats: 

1 dyn ~ E ■ K' 1 ) 

Pi 
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VI. THE THERMODYNAMIC QUANTITIES FOR THE COUPLED SPIN FLUID 

• The heat absorbed or delivered by the heat engine 
Using Eq.([T5|) the heat Qh/c absorbed or delivered becomes: 

Qi = (exp(-rr i )-l)(a; i b 1 + Jb 2 ) (72) 

where i = h/c 

• The work absorbed or delivered by the heat engine 

The power, Eq. (flljl. is (^) = Bx(t)u. Therefore, the work becomes: 

W = biwdt . (73) 

J ri 

• Entropy production. 

The entropy production per cycle, VS cyc [ e , created on the boundaries becomes (Cf. 
Eq. jnj): 

VS cycle = - (Q AB /T h + Qcd/T c ) . (74) 

• Efficiency. 

The efficiency per cycle, r] cyc[e is: 



f Tf bicjdt , 

Vcycle = W Q AB = - J " \ . 75 

(exp(-rri) - l)(^ibi + Jb 2 ) 
The maximal efficiency of the engine is: 



- i v 7 ^ 7 ^ 

The upper bound should be the Carnot's efficiency, a bound correct for all J and the 
fact that ^ > 2k : 

, ,2 I 72 , ,2 T 2 

1 - ^ + J < l-^L < i_l£. (76) 
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VII. THE CYCLE OF OPERATION: THE OTTO CYCLE 



The operation of the heat engine is determined by the properties of the working medium 
and by the hot and cold baths. These properties are summarized by the generator of the 
dynamics C The cycle of operation is defined by the external controls which include the 
variation in time of the field with the periodic property u(t) = u(t + r) where r is the total 
cycle time synchronized with the contact times of the working medium with the hot and cold 
baths Th and r c . In this study a specific operating cycle composed of two branches termed 
isochores where the field is kept constant and the working medium is in contact with the 
hot/cold baths. In addition two branches termed adiabats where the field u(t) varies and 
the working medium is disconnected from the baths. This cycle is a quantum analogue of 
the Otto cycle. 

The dynamics of the working medium has been described in Sec. IIHI The parameters 
defining the cycle are: 

• Th and T&, the hot/cold bath temperatures. 

• I\ and r c , the hot/cold bath heat conductance parameters. 

• 7h and 7 C , the hot/cold bath dephasing parameters. 

• J-the strength of the internal coupling 

The external control parameter define the four strokes of the cycle (Cf. Fig. [2J): 

1. Isochore A —>■ B: when the field is maintained constant, uj = Uf,, the working medium 
is in contact with the hot bath for a period of r^. 

2. Adiabat B — > C: when the field changes linearly from Ub to u a in a time period of Tb a - 

3. Isochore C — > D: when the field is maintained constant u = u a the working medium 
is in contact with the cold bath for a period of r c . 

4. Adiabat C — > A: when the field changes linearly from u a to uj^ in a time period of t &. 

The trajectory of the cycle in the field and the entropy plane (u,Se) is shown in Fig. |21 
employing a numerical propagation with a linear uj dependence on time. 
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FIG. 2: The heat engine's optimal cycles in the (uj,Ss) plane. The upper red line indicates the 
energy entropy of the working medium in equilibrium with the hot bath at temperature T^ for 
different values of the field. The blue line below indicates the energy entropy in equilibrium with 
the cold bath at temperature T c . The cycle in green has an infinite time allocation on all branches. 
It reaches the equilibrium point with the hot bath (point E) and equilibrium point with the cold 
bath (point F). The inner cycle ABCD is the optimal cycle with the optimal time allocation on all 
branches, calculated numerically for a linear ui dependence on time, t/j = 3.0108 r^ a = 0.301, r c = 
3.014 r c h = 0.346. The external parameters are: uj c = 5.382, uj^ = 12.717, J = 2., Th = 7.5, T c = 
1.5, F h = 0.382, T c = 0.342, j h = 7c = 

A different perspective on the dynamics during the cycle of operation is shown in Fig. 
EJ displaying the cycle trajectory in the bi,b 2 ,b 3 coordinates. The hypothetical cycle with 
infinitely long time on all branches would include the equilibrium points E and F. The cycle 
trajectory is planar on the B 3 = plane as can be seen in panel C. The cycle ABCD 
with finite time allocation spirals around the infinitely long time cycle with an incursion 
into the B 3 directions. The reference cycle with infinite time allocation on all branches is 
characterized by a diagonal state p e in the instantaneous energy representation. The slow 
motion on the adiabats allows the state p to adopt to the changes in time of the Hamiltonian, 
which therefore can be termed adiabatic following. If the time allocation on the adiabats is 
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FIG. 4: Three cycles of operation based on the analytic solution in the (uj,Se) plane. The orange 
inner cycle has the shortest time allocations (jh = 2. T ba = r a {, = 0.05, r c = 2.1). The green 
cycle shows the corresponding (lo,Svn) plot. The magenta cycle has longer time allocations r/j = 
r c = 15. Tb a = T a t, = 0.015, while the black cycle has infinite time allocations on all branches 
therefore Se = Svn- This cycle touches the isothermal equilibrium points E and F. The common 
parameters for all the cycles are: J = 2.,r = 0.96, Th = 7.5, T c = 1.5, = r c = 0.3243, -jh = 
7c = o,uj a = 5.08364, oo b = 11.8675. 

short, non-adiabatic effects take place. In the sudden limit of infinite short time allocation 
on the adiabat, the state of the system has no time to evolve p(t{ + r a b) = p(tj). The 
Hamiltonian will then change from Hi = (^(t^Bi + JB 2 to Hf = u(ti + r a i ) )'B 1 + JB 2 
therefore the representation of the state p e {t% + t^) in the new energy representation is 
rotated by an angle 9 = (8i — Of) compared to the former one. Where, 9i = arcsin( J/Q(ti)) 
and 9f = arcsin ( J/Q(ti + r a b)). When following the direction of the cycle, the energy- 
entropy increases on the adiabts. This is evident in both Fig. Eland Fig. |3J This entropy 
increase is the signature of nonadiabatic effects reflecting the inability of the population on 
the energy states to follow the change in time of the Hamiltonain. As a result the energy 
dispersion increases. Since the evolution on these branches is unitary, Svn is constant. 
When more time is allocated to the adiabats the increase in Se is smaller. For infinite time 
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allocation Se = Svn- In this case the state of the woking medium is always diagonal 
in the energy representation. The larger curvature of the entropy increase in the analytic 
result of Fig. compared with the numerical result of Fig. |21 reflects the difference in 
the dependence of u(t) on time. When the analytic functional form of u(t) is used in the 
numerical propagation the numerical solution converges to the values of the analytic solution. 
This convergence test was used as a consistancy check for both methods. Convergence was 
not uniform for all elements in the propagator (Cf. Eq. PU|) and Eq. ()46|) ). Comparing 
the elements of the numerical propagator U a (r ab ) to the elements of analytic U a (r ab ), showed 
that the largest discrepency between the individual elements at t = T ab was less than 1CT 3 
when a time step of At = r ab /1000 was used. 

In Fig. |S]the cycle of operation is presented in the energy-entropy internal-temperature 
coordinates (SE,Td yn )- The cycles shown corresponds to the analyticial cycles of Fig. HI 
The discontinuities in the short time cycle reflect over-heating in the compression stage as 
shown as the difference between the point A and A' in Fig. The heat accumulated is 
quenched when the working medium is put in contact with the hot bath. This phenomena 
has been identified in measurements of working fluid temperatures in actual heat engines or 
heat pumps [26j. A discontinuity as a result of insufficient cooling of the woking medium 
in the expansion branch is also evident in the short time cycle. The magnitude of these 
discontinuities is reduced at longer times and dissapear for the infinite long cycle where the 
working fluid reaches thermal equilibrium with the hot bath at point E and with the cold 
bath at point F. In this case both adiabatic branches are isoentropic. It is clear from Fig. 
that for the cycles with vertical adiabats the work is the area enclosed by the cycle trajectory. 
When the time allocation on the adiabats is restricted this is no longer the case since due to 
the entropy increase, the area under the hot isochore does not cover the area under the cold 
isochore. Additional cooling is then required to dissipate the extra work required to drive 
the system on the adiabats at finite time. 

VIII. THE EFFECT OF PHASE AND DEPHASING. 

The performence of the heat engine explicitly depend on heat and work which constitute 
the energy Do other observables, incompatiable with the energy, influence the engins 
performence? Examining the cycle trajectory on the isochores in Fig. EJ in addition to the 
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FIG. 5: The cycles in (Sg, T^) planes. The inner cycle A,B,C,D corresponds to the short time 
cycle of Fig. (jljl. The magneta cycle is the long time cycle and the black cycle H,E,G,F corrosponds 
to the cycle with infinite time allocation on all branches. The rectangle, including points I,E,K,F 
is the work obtained in a Carnot cycle operating between Th and T c . The shaded area H,E,G,F 
represents the maximum work of the Otto cycle. The area below the A,B segment is the heat 
trasfered from the hot bath Q^. The area below the D,C segment is the heat transfered to the cold 
bath Q c . 

motion in the energy direction, towoard equilibration, spiraling motion exists. This motion 
is characterized by amplitude and phase of an observable in the plane perpendicular to the 
energy direction. The phase of this motion advances in time, i.e. <fi oc t. The concept 
of phase has its origins in classical mechanics where a canonincal trnsformation leads to a 
new set of action angle variables. The conjugate variable to the Hamiltonain is the phase. 
In quantum mechanics the phase observable has been a subject of contineous debate 
For a harmonic oscillator it is related to the creation and anhilation operator a 39, 4(|. In 



analogy the raising/lowering operator is denned: 

J_(_jB 1+ cuB 2 ±^B 3 ) , (77) 



V2n 
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which has the following commutation relation with the Hamiltonian: 

[H, L±] = ±y/2QL± . (78) 

The free evolution of L + therefore becomes: L+(t) = e lv ^ f2t L + (0) which defines the phase 
variable through: (L+) = re^, therefore <p = arctan f z^+^ V ^ corroboration for this 
interpretation is found by examining the state p e in the energy representation (Cf. Eq. (|63j). 
The off diagonal elements are completely specified by the expectation values of L±. 

The dynamics of L± on the isochores includes also dissipative contributions which can 
be evaluated using Eq. (|3^j): 

L± = ± iV2ttl± - (r + 2 7 fi 2 ) L± (79) 

Examining Eq. (|79j) it is clear that the amplitude of L± decays exponentially with the 
rate ^- = T + 2-ffl 2 , where T is the dephasing contribution due to energy relaxation and 
t^t = 2 r yQ 2 is the pure dephasing contribution. 

Both Fig. El and Fig. |H1 show that the dephasing is not complete at the end of the 
isochores. A small change in the time allocation in the order of 1/Q can completely change 
the final phase on the isochore and on the initial phase for the adiabat. This means that the 
cycle performance characteristic becomes very sensitive to small changes in time allocation 
on the isochores. This effect can be observed in Fig. [7| for the power and Fig. |H] for the 
entropy production. Examining Fig. [7|reveals that increasing J increases the "phase" effect. 
For J = 2 for specific time allocations the power can even become negative. Increasing the 
dephasing rate either by adding pure dephasing or by changing the heat transfer rate reduces 
the "noise". This can also be seen in Fig. |S| An interesting phase effect can be observed 
in Fig. El where the cycle is displayed in the (SE,Td yn ) plane. The inner (solid black) cycle 
shows an energy-entropy decrease in the compression adiabat. The reason for this decrease 
is a phase memory from the compression adiabat which is due to insufficient dephasing on 
the cold isochore. Additional pure dephasing eliminates this entropy decrease as can be seen 
in the dashed black cycle. This cycle is also pushed to larger entropy values. The orange 
cycles are characterized by a longer time allocation on the isochores. For these cycles the 
energy-entropy always increases on the adiabats. This cycle is shifted by dephasing to lower 
energy-entropy values. 
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FIG. 6: The modulus and phase of h± as a function of time. The dashed lines include additional 
pure dephasing (7^ = 0.01, 7 C = 0.03). The common parameters are: T^=7.5, T c = 1.5,1^ = 
r c = 0.34, uif, = 11.8675, w a = 5.083, The total cycle time is r = 2.4 where, t/j = r c = 1, 
r ba = 0.2,T ab = 0.2. 

IX. DISCUSSION 

Quantum thermodynamics is the study of thermodynamical phenomena based on quan- 



tum mechanical principles 43]. To meet this challenge, quantum expectation values have 



to be related to thermodynamical variables. The Otto cycle is an ab-initio quantum model 
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FIG. 7: The power produced by the engine as a function of the time allocation on the hot isochore. 
For the green cycle J=l and = r c = 0.324. For the blue cycle J=2 and = r c = 0.324. 
For the red cycle J=2 and Th = T c = 0.162. The three colored cycles have no pure dephasing 
lh = 7c = 0. With addition of dephasing 7^ = 0.01 and 7 C = 0.03 the "noise" is eliminated and 
the three cycles collapse to the solid black lines. The common parameters are: 7^=7.5, T c = 1.5, 
Ub = 12.717, w a = 5.382, The total cycle time r is:=6.74, T& a = 0.3, T a b = 0.34. 

for which analytic solutions have been obtained. The principle thermodynamical variables: 
energy entropy and temperature are derived from first principles. The solution of the quan- 
tum equations of motion for the state p, enables tracing the thermodynamical variables for 



each point on the cycle trajectory. T 



finite-time-thermodynamics 



lis dynamical picture supplies a rigorous formalism for 
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An underlying principle of finite-time-thermodynamics is that operation irreversibilities 
are inevitable if a process is run at finite rate. Moreover these irreversibilities are the source 
of performance limitations imposed on the process. The present Otto cycle heat engine in 
line with FTT is subject to two major performance limitations: 

• Finite rate of heat transfer from the hot bath to the working medium and from the 
working medium to the cold bath. 
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FIG. 8: Entropy productions T>S cyc i e Eq. (|74|) . as a function of the time allocation on the hot 
isochore. The notations are the same as Fig. [7| 

• Additional work invested in the expansion and compression branches is required to 
drive the adiabats at a finite time. 

The finite rate of heat transfer limits the maximum obtainable power V The present 
Otto engine model is not an exception, showing similarities with previous studies of discrete 
quantum heat engines jillfi, 10 . fll| . 

The irreversibility caused by the finite time duration on the adiabats is the novel finding 
of the present study as well as the preceding short letter (^J. This irreversibility is closely 
linked to the quantum adiabatic condition. The nonadiabatic irreversibility is caused by 
the interplay of the noncommutability of the Hamiltonian at different points along the cycle 
trajectory and the dephasing caused by coupling to the heat baths on the isochores. In 
the present Otto cycle these contributions are separated in time. The non-adiabaticity can 
be characterized by an increase in the modulus of (L±) on the adiabats. Dephasing, i.e. 
exponential decay of the modulus of (L±) is induced by the coupling to the baths on the 
isochores. 

The dynamics of the Ey. operator associated with the phase can be compared to the B± 
operator associated with the internal correlation between the spins (Cf. (JB2J)). The absolute 
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FIG. 9: The influence of dephasing on the cycle of operation in the {SE,Tdy n ) plane. Solid curves 
correspond to an operation without pure dephasing, The dashed curves represent cycles including 
pure dephasing. For the black cycles the time allocations on the isochores are: Th = t c = 0.6. 
The pure dephasing parameter is jh = 7c = for the solid lines, and 7^ = 0.005, j c = 0.015 for 
the dashed lines. For the red cycles the allocated times on the isochores are: = 2.,r c = 2.1 
with 7/j = 7 C = for the solid lines, and 7^ = 0.01, 7 C = 0.03 for the dashed lines. The common 
parameters for all four cycles are: J = 2.,Th = 7.5,T C = 1.5,^ = r c = 0.3243, r a b = Tba = 0.015. 

value of |B±| oscillates on all branches of the cycle never reaching zero. This is not surprising 
since B-t does not commute with the Hamiltonian. The "angle" <$b = arctan(&3/62) is excited 
for small cycle times. For cycles with large time allocation on the isochores, 4>b is found to 
be close to zero. These observations reflect the two types of correlations between particles. 
A "classical" correlation and a quantum correlation meaning EPR 0, 4^ entanglement 
between particles. The general trend is therefore for the engine to become more "classical" 
when the cycle times become longer. In this case the state follows the energy direction and 
in addition entanglement between particles is small. Adding pure dephasing has a similar 
effect. A continuous measurement of energy during operation will also lead to effective pure 
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dephasing. For short cycle times quantum effects become important. The entropy decrease 
on the adiabats which is the result of "phase" memory is such an example. The quantum 
effect which influences the performance is the excess work on the adiabat due to the inability 
of the state to follow the energy direction. 
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APPENDIX A: THE F OPERATORS 



The method of construction of Fj is based on identifying the operators with the raising and 
lowering operators in the energy frame. The matrix C which diagonalizes the Hamiltonian 
becomes: 
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iL n n J- \ 

V2 u u V2 













\ji o o -fj 



( \ 






\ 7SJ 



(A2) 



f -fi x \ 
10 
10 

\ x o uj 

The down transition rates k j are chosen to be equal for all the four transitions, while the 
raising transitions k j comply with detailed balance. Schematically the eight transitions are: 

FT? T? T? T? T? T? T? 
1 i 2 *3 *4 *5 *6 *7 *8 

E/2 E2 E3 E3 E4 E4 E4 E4 

Ei E\ Ei Ei E2 E2 E% E3 



(A3) 
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a. Detailed presentation of a few Fi operators. 
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The F operator for the transition E\ to E 2 , is Fi^+2 = Fi; In the energy picture, it is 
simply: 
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Using the matrix C to transform back to the polarization picture leads to: 
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Using a similar procedure all the Fi in the polarization picture become: 



Fx = Fx. 
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